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SUPERSONIC LIFT AND PITCHING MOMENT OF THIN SWEPTBACK 
TAPERED WINGS PRODUCED BY CONSTANT 
VERTICAL ACCELERATION 

SUBSONIC LEADING EDGES AND SUPERSONIC TRAILING EDGES 
By Frank S. Malvestuto, Jr., and Dorothy M. Hoover 

SUMMARY 


Theoretical approximations to the nondimensional lift and pitching 
moment produced by constant vertical acceleration, the Cl^ and C^ 

derivatives, respectively, are derived for a series of thin sweptback 
tapered wings with streamwise tips. The analysis is essentially the 
application of a recently published solution of the linearized time- 
dependent wave equation for wings in accelerated motion. The results 
are independent of camber and thickness and are applicable for a r ang e 
of supersonic speed for which the wing is wholly contained between the 
Mach cones springing from the wing apex and. from the trailing edge of 
the root section (subsonic leading edge and supersonic trailing edge). 

Design curves are presented which permit rapid estimation of the 
derivatives Cl^ and Cm^ for given values of aspect ratio, taper 

ratio, Mach number, and leading-edge sweep. 


INTRODUCTION 


The formulations of linearized supersonic aerodynamics ha ve allowed 
the theoretical derivation of many of the important longitudinal- and 
lateral-stability derivatives of various wing configurations. Recently, 
attention has been focused on the sweptback tapered wing with wing' tips 
parallel to the wing plane of symmetry (hereinafter referred to as the 
"sweptback tapered wing"). Available stability derivatives for this 
wing for a wide range of Mach number now include the lift-curve slope Ci^, 

references 1 to 3; the damping-in-roll derivative C j p , references 2 to 4; 

the lateral-force and yawing-moment derivatives Cy^ and C n , 

P P 
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respectively, references 5 and 6 ; and the longitudinal derivatives Cn^ 
(pitching moment due to angle of attack), CL q (lift due to pitching), 
and C mq (damping in pitch), reference 7 . 

The present paper is an extension of the previous investigations 
dealing with the sweptback tapered wing traveling at supersonic speeds 
and is concerned with the lift and the pitching moment resulting from 
constant vertical acceleration, the Cj^ and derivatives, respec- 

tively. Theoretical expressions are derived herein which approximate 
with practical accuracy these -two acceleration derivatives. 

Basically, the analysis depends upon, the solution of the linearized 
time -dependent wave equation for wings in accelerated motion reported by 
Gardner and by Watkins in references 8 and 9 , respectively, and the appli- 
cation thereof to triangular wings in reference 10. The analyses of refer- 
ences 8 and 9 have in effect demonstrated that for a wing the time- 
dependent potential for constant vertical acceleration may be compounded 
of two time-free or steady-state potentials, one for angle of attack and 
the other for steady-state pitching. In terms of forces and -moments, 
the analysis of reference 10 shows that such a decomposition of the time- 
dependent potential allows the derivatives and to be 

expressed in terms of the steady-state -pitching and angle -of -attack 
derivatives (CL q , and C m a ) P lus additional terms dependent 

upon the surface velocity potential and corresponding pressure function 
for constant angle of attack. The results of reference 10 and the 
present results for the Cp^ and derivatives are restricted to 

relatively slow rates of acceleration such as arise in the study of air- 
plane stability. Approximate expressions for the derivatives Cj^, Cm^ 

Cp q , and. C mq for use in the evaluation of the derivatives of sweptback 

wings contained in this paper have previously been determined in refer- 
ences 3 and 7* Approximate expressions for the terms dependent upon the 
surface velocity potential and corresponding pressure distribution are 
derived herein. 

. The results of the analysis are applicable for a range of Mach 
number which allows the leading edge to be subsonic and the trailing 
edge to be supersonic with the additional trivial limitation that the 
Mach 'lines from the leading edge of the wing tip cannot intersect on 
the wing or intersect the opposite wing edges. Design curves are pre- 
sented which permit rapid estimations of the derivatives Cp^ and Cm^ 

for given values of aspect ratio, taper ratio, Mach number, and leading- 
edge sweep . 
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SYMBOLS 


x,y,z Cartesian coordinates of an arbitrary point 

u i> v i> w i induced flow velocities along x, y, and z body axes, respec- 
■ tively (see fig. 2(a)) 


u,v,v 

X,Y,Z 

V 

d 

a 

<1 

M 

B 


incremental flight velocities along x, y, and z stability 
axes, respectively (see fig. 2(b)) 

forces parallel to x-, y-, z-axes, respectively 

flight speed 

angle of attack (w/v) 

rate of change of a with time (da/dt) 

pitching velocity about y-axis 

stream Mach number (V/Speed of sound) 

Mach angle 

cotangent of Mach angle 


/M 2 - 1 


X 

A 

9 0 = tan e 
5 

tan e 
tan (i 


angle between leading edge and axis of wing symmetry (see 
fig. 1(a)) 

taper ratio of wing 


m = 


leading-edge sweep (90° - e) 

angle between trailing edge and axis of symmetry (see fig. 1(a)) 
= 0 n B 


0) 


geometric parameter of wing 


'2c r 0 o 


4e o km 

A(1 + X) AB(l + XT j 


n = Hirl = 1 - (i - 


= _x_ 

X0Q 


V 
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x 

c r 

•c 


S 

A • 

* 

X 

0 

t 

a 

P 

AP 

<*V 

C P 

M' 


wing span 

x-component of center -of -gravity position 
wing root chord 
over -all length 


mean aerodynamic chord l — 


nb/2 


(Local chord) dy = 


2c r (la> 2 - 3a)(l - n) + (1 - n) 2 ^]\ 
3 cq 2 (1-+ \) / 


3a> (l + X) 

distance of leading edge of local wing chord behind apex of 
wing 

wing area 

aspect ratio of wing 

steady-state potential corresponding to a unit pitching 
velocity about y-axis 

steady-state potential corresponding to unit angle of attack 
perturbation velocity potential on upper surface of wing 
time 

speed of sound 
density of air 

local pressure difference between lower and upper surfaces of 
airfoil, positive in sense of lift 

lift distribution for unit pitching velocity 
lift distribution for unit angle of attack 


/ 


pressure coefficient 
pitching moment 


AP 
1 „2 
£P V / 
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E'(m) complete- elliptic integral of the second kind with modulus k 



When x, y, and t are used as subscripts, the respective partial 
derivative is indicated. For example. 



Primed symbols X', AP', and 0' refer to the wing-tip region shown 
in figure 1. 


ANALYSIS 

Scope 


The sweptback wings considered in this paper are sketched in fig- 
ure 1. In the following analysis, the plan form with sweptback trailing 
edge (fig. 1(a)) is generally considered to be the typical wing, but the 
results of the analysis are equally valid for the wing plan form with 
the sweptforward trailing edge (fig. 1(b)). The orientation of the wing 
with respect to a body system of coordinate axes used in the analysis is 
indicated in figure 2(a). The surface velocity potential, the basic 
pressure distribution, and the stability derivatives are derived with 
respect to this system. Figure 2(b) shows the wing orientated with 
respect to the stability_axes system with the origin of the system at 
the arbitrary location (x, 0, 0) rearward of the apex of the wing. The 
transformation formulas that allow the determination of the derivatives 
with respect to the stability system once they are evaluated with 
respect to the body axes system are presented in table I. 

The analysis is limited to wings having small thickness and camber 
and that are not yawed with respect to the free-stream direction. The 
derivatives are valid only for a range of supersonic speed which allows 
the leading edge to be subsonic and the trailing edge to be supersonic. 
The terms "subsonic leading edge" and "supersonic trailing edge" refer 
to the conditions that the Mach number of the stream component normal to 
the leading edge is less than 1 and that the Mach number of the stream 
component normal to the trailing edge is greater than 1, respectively. 

An additional trivial restriction is that the Mach lines emanating from 
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the wing tips cannot intersect on the wing. These conditions, expressed 
mathematically as restrictions for B cot A, are as follows: For 

BA(l + A.) = 2 

BA(l + X.) . „ . . ^ , 

BA( 1 + X) + 4(1 - X) = B cot A = : 

,and for BA(1 + X) < 2 


BA( 1 + X) < _ . . . 

BA(1 + X) + 4(1 - X) = B COt A 


< BA(1 + X) 

" 4 - BA(1 + ,X) 


Basic Considerations 

The derivation of the derivatives Cl^ and Basically depends 

upon the solution of the linearized potential equation for time- 
dependent motion 


B ^xx “ ^yy + . 2 ^xt + 2 ^tt " 0 (1) 

a ' a 

required to fulfill the boundary condition on the wing (approximately 
in the z = 0 'plane) 



It was reported in the section of the analysis of reference 10 
dealing with the acceleration derivatives of triangular wings that 
Gardner has, in effect, shown that a proper solution of equation (l) is 




a 


B 


2 . 


*•- 



( 2 ) 


where is the steady-state potential corresponding to a unit pitching, 
velocity about the y-axis and X is the steady-state potential corre- 
sponding to a unit angle of attack. Thus, the solution by Gardner and 
an essentially equivalent solution by Watkins (reference 9) allow the 
time -dependent potential for an angle of attack at to be expressed in 
terms of two time -independent or steady-state potentials, one for a 
constant angle of attack and the other for steady pitching. The appli- 
cation of the solution expressed in equation (2) permits with relative 
ease the evaluation of the pressure distribution for vertical motion 
with constant acceleration. The corresponding derivatives CL& and Cm^ 
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may then be derived rather simply in terms of integral expressions of 
the surface velocity potential and surface pressure distribution for 
constant unit angle of attack together with terms involving the non- 
dimensional forces and moments associated with angle of attack and 
steady pitching. The general formal ievelopment of the derivatives Cl^ 

and Cn£ based upon the velocity potential expression of equation (2) 

has been carried out in reference 10. Adhering to the analysis in that 
paper, the lifting-pressure distribution at time t = 0 for the angle * 
of attack at is obtained from the surface velocity potential by 

AP = 2p(v0 x . + 0 t ) t=Q 


.= 2pVd 




VB^ 



B 


a 

2 


M 2 x 


^(AP)„ =1 - — (AP)^ - 2pX 


'q=l 


( 3 ) 


where (AP)^_^ is the lift distribution for' unit pitching velocity 
about the y-axis and (AP) a _^ is the lift distribution for unit angle of 

attack. The choice of time t = 0 eliminates the lift due to angle of 
attack and leaves only the increment due to time rate of change of angle 
of attack. 


Integration of equation ( 3 ) to obtain the lift and moment and 
reduction to coefficient forms yields 



where the integrations are performed over the wing plan form. For 
convenience, the following symbols will be used: 

- S + 2C -a 
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N 2 'iiTl ds 

so that 

CL a = ^( Cl <*)i “ ^( Cl o ) 2 


and, similarly, 



so that 



Satisfactory approximate expressions for the derivatives C m a , CLq, 
and- that occur in equations (4) and (5) for Cl^ and 

respectively, have previously teen determined and presented in refer- 
ence 7 tor the wing plan forms and Mach number range considered in this 
paper. The double integral terms remain to be evaluated for a complete 
estimation of C]^ and C^. These integrals and the approximate 

expressions therein for the potential X and the pressure (AP) a= ]_ are 
considered in the following sections. 


Expressions for the integrals 



and 


X y dS « 


An approxi- 


mation to -the steady-state surface velocity potential of a sweptback wing 
at a unit angle of attack has been derived in reference 3- The entire 
approximation occurs in the expression for the surface velocity potential 
for the region of the wing internal to the Mach cones from the leading 
edge of the wing tips. For the region of the wing external to the wing- 
tip Mach .cones, 'the velocity potential distribution is exact and is the 
same as the linearized potential distribution for the corresponding part 
of a triangular wing for the same leading-edge sweep and Mach number. 
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The surface velocity potential for the triangular wing may he derived 
from the analysis of reference 11 and is expressed as follows for the 
upper surface of the wing for unit angle of attack: 


v 0 V\|e 0 2 * 2 - y 2 

X = a E'Tm) 


or 


X = 


V0 o x\|l - 
E ' (m) 


( 6 ) 


where V is a conical coordinate y/0 Q x proportional to the slope of 
a ray from the apex of the wing and E'(m) is the complete elliptic 

integral of the second kind dependent upon m = *'. an e . The variation 

tan n 

of with m is presented in figure 3- From reference 3, the 

corresponding approximate expression for the potential of the wing-tip 
region is given by 


Y , = £_ = 2V / 2(mx + By)(b - 2y) 
a tc ]j B(l + m) 


(T) 


The approximate surface velocity potential for the entire wing, for a 
unit angle of attack is then the sum of the expressions (6) and ( 7 ). 
The integral terms of equations (4) and (5) containing X may now be 
expressed in the following form (see fig. 4 for regions): 


rr x ^ , i6e ° 

B%c JJ wing V B^ScE ' (m) 


If X] F~ 


v ds + 


Region 

oegh 


32\/2~ 


nB 2 Sc^B(l + m) 


\f( mx + By)(b 
Region Y 
efg 


- 2y) dS (8) 
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and 


8 


1 . 


B 2 Sc 2 00 Wing 


x | dS = - 

v •oCQ-^l? 1 I 


B^Sc^E ' (m) 


x 2 \/l 


Region 

oegh 


\pT- 


p 

v* as + 


32 \[ 2 ~ 


itB 2 Sc 2 ^B(l. + m) 


Region 
efg 


x\/(mx + By)(b - 2y) dS 


(9) 

The quasi -conical nature of the integrands (that is, of the form x n f(v) 
of the first integrals of the right-hand side of equations (8) and (9) 
allows these integrals to be easily evaluated by use of the "triangular" 
integration procedure considered in reference 12. When this procedure 

is used, the integrand xJl - v 2 dS may be written as 


Region ohgo 

i c 3 e £ ZZ 

3 ' ° (1 - nv) 3 


dv 


Region ogeo 
1 (l + m)^b^ \/l - y 2 


2K 


(l + mv)^ 


dv 


and, similarly, the integrand x 2 \/l 
Region ohgo 


v 2 dS becomes 

Region ogeo 


1 

5 


c r 



\/l - V 2 

(1 - nv)^ 


dv 


1 

64 


(l 


+ m)HA ' Vl - v 2 
9 q 3 (l + mv) 1 *' 


dv 
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Expression for the integral 



dS.- The evaluation of 


a=l 


the integral 



dS necessitates a knowledge. of the lifting 


/ a=l 

pressure for the wing at an angle of attack. For the region of the wing 
external to the wing-tip Mach cones, the pressure distribution is merely 
the exact linearized pressure distribution for the corresponding part of 
the triangular wing for the same value of leading-edge sweep and Mach 
number. From reference 11 , the following expression may be obtained for 
the lifting-pressure coefficient for a unit angle of attack 



40 o 2 x 


E'(m)\/e 0 2 x 2 - y 2 


or, in terms of the conical coordinate v = ^ , 

0Q X 


■ AP 
1 ,„ 2 
V2 


40 , 


pV / a=1 E'(m)/l - v 2 


(14) 


The corresponding expression for the approximate lifting-pressure distri- 
bution for the region of the wing internal to the wing-tip Mach cones has 
been derived in reference 3 . From therein, the lifting-pressure coeffi- 
cient for unit angle of attack is as follows: 



Substitution of the 
tions (14) and (15) 


sum of the component pressures expressed by equa- 
' AP \ 

of the integral to be evaluated yields 

a=l 


for 
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RESULTS AND DISCUSSION 


The preceding section involved the evaluation of the acceleration 
derivatives Cl^ and as functions of the steady-state pitching 

derivatives C^, C m ^j determined in reference 7. and the several 

integral terms determined herein which are functions of the surface 
velocity potential and the corresponding pressure distribution for con- 
stant unit angle of attack. The evaluation of these derivatives are . 
approximate in that the contribution of the wing-tip region was obtained ' 
by use of an approximate surface velocity potential derived in refer- 
ences 3 and 7 by use of the point -source method of Evvard. 

An indication of the lifting-pressure distribution over a typical 
wing (sketched in fig. 5) for three Mach line configurations can be 
obtained from figure 6 . Figure 6 (a) shows the pressure distribution for 
the case of a subsonic leading edge and a supersonic trailing edge; 
figure 6 (b) , for a subsonic leading edge and a sonic trailing edge; and 
figure 6 (c), for a sonic leading edge and a supersonic trailing edge. 

As expected, a pronounced finite drop in the resultant lifting pressure 
occurs across the inboard Mach line from the wing tip. The resultant 
pressure distribution is generally negative for the wing plan form con- 
sidered. For the Mach line configuration for which the leading edge is 
sonic and the. trailing edge is supersonic (fig. 6 (c)), the pressure along- 
section a-a becomes slightly positive in the wing-tip region, in particu- 
lar in the vicinity of the trailing edge. Results of additional computa- 
tions for the variation of pressure loadings with aspect ratio indicate 
that, as the aspect ratio is decreased, the resultant pressures tend to 
become positive and eventually give rise to a positive (lifting) force 
(C 14 ) and a stable nose-down pitching moment whereas, for the 

wing considered for figure 6 , it is obvious that the integrated pres- 
sures will produce a negative Cp^ and an unstable C^. 

The exact linearized solution for the lifting pressure in the wing- 
tip region induced by constant vertical acceleration is not available at 
present although extremely accurate approximations to the exact linearized 
solution that require laborious calculations may be obtained by the 
methods of reference 13. The present approximate method should be ade- 
quate for determining the integrated lift and pitching moment, especially 
since the pressures are relatively very small in the tip' region where 
the approximations of the method apply. For a practical evaluation of 
the derivatives Cp^ and it is believed that the wing-tip regions 

may be completely neglected. 

A series of generalized curves that allow rapid estimation of the 
derivatives Cp^ and Cm^ is presented in figures 7 and 8 , respectively. 



18 


NACA TN 2315 


For specified values of aspect ratio, taper ratio, Mach number, and 

leading-edge sweep, the derivative Ct. = ~ [Ct.^ - — \ is esti- 

^ B 2 V a /l B 2 \ a /2 

mated from figure 7 and the derivative ~ 

estimated from figure 8. In figures 7 and 8, the segmental dashed parts 
of the design curves to the left of the boundary lines labeled 
"Sonic T.E." should be disregarded for estimations; they correspond to 
conditions not treated in this paper, for which the trailing edge is 
subsonic (Mach lines ahead of trailing edge) and hence allows trailing - 
edge disturbances to affect the part of the wing bounded by the Mach 
lines from the trailing-edge apex and trailing edge itself. These 
dashed parts have been presented, however, to indicate the trend of the 
variations and to act as an upper limit below which the true values of 
the derivatives would lie for configurations for which the trailing edge 
is subsonic. The boundary lines labeled "Sonic L.E. " indicate the 
other limit of validity for which the Mach lines coincide with the 
leading edge . 

It should be noted that the derivatives determined from the design 
curves of figures 7 and 8 are with respect to a set of axes located at 
the apex of the wing. The derivatives with respect to an arbitrary 
center -of -gravity location (x = x, y = 0, z = 0) in the stability axes 
system may be easily obtained by transformation formulas presented in 
table I. 

Specific variations of the derivatives Cl^ and (in the sta- 

bility axes system) with each of the parameters - aspect ratio, taper 
ratio, Mach number, and leading-edge sweep - are presented in figures 9 
and 10 , respectively. The variation of Cn^ with aspect ratio plotted 

in figure 10 shows the interesting point, mentioned previously in the 
discussion of the pressure distributions of figure 6, that the pitching 
moment which is generally unstable at the higher aspect ratios becomes 
stable as the aspect ratio is decreased (below aspect ratio of approx. 2.25 
for the illustrative case shown) . Similar results have been noted for 
the rectangular wing (reference 14 ), the triangular wing (reference 10), 
and the sweptback wing of zero taper ratio (reference 12). 


CONCLUDING REMARKS 


Theoretical approximations to the nondimensional lift and pitching 
moment produced by constant vertical acceleration, the Cl^ and Cm^ 

derivatives, respectively, have been derived for a series of thin 
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sweptback tapered wings with streamwise tips. Results are applicable 
for a range of supersonic speed for which the wing is wholly contained 
be tween the Mach cones springing from the wing apex and from the trailing 
edge of the root section. 

The results presented herein for the derivatives Cl^ and C ^ are 

applicable within the limitations of the linearized theory since the 
approximation to the exact linearized solution in the wing-tip region 
is satisfactory and this region has a minor effect on the values of the 
derivatives for the entire wing. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 22, 1950 
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Figure k.- Sketch of right panel of wing and associated data defining 
limits of integration for the evaluation of forces and moments. 
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Figure 7.- Variation of Bfi^ and B^fc^ with B cot A for various 

values of AB and X. Body axes system, origin at apex of wing. 

Dashed portions of curves have limited significance (see "Results and 
Discussion" ) . 
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Figure 7* - Continued, 
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Figure 7-- Continued 
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Figure 7 .- Continued. 
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(c) X = 0.75. 


Figure 7 .- Continued 
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(c) Concluded. A. = 0.75* 


Figure 7 .- Concluded, 







values of AB and X.. Body axes system, origin at apex of wing. 
Dashed portions of curves have limited significance (see "Results 
Discussion" ) . 
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Figure 8.- Continued 
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(c) X = 0.75. 
Figure 8.- 


Continued. 
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(c) Concluded. X. = 0.75* 
Figure 8.- Concluded. 












